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Abstract 



The role of higher Fock-state contributions to exclusive charmonium decays will be dis- 
cussed. It will be argued that the J/ip (VO decays into baryon-antibaryon pairs are dominated 
by the valence Fock-state contributions. P-wave charmonium decays, on the other hand, re- 
C ■ ceive strong contributions from the ccg Fock states since the valence Fock-state contributions 

are suppressed in these reactions. Numerical results for the decay widths of J/ip(ip') — > BB 
■ and Xcj — ► PP will be also presented and compared to data. 
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1 Introduction 



Exclusive charmonium decays have been investigated within perturbative QCD by many au- 
thors, e.g. Q. It has been argued that the dominant dynamical mechanism is cc annihila- 
tion into the minimal number of gluons allowed by colour conservation and charge conjuga- 
tion, and subsequent creation of light quark-antiquark pairs forming the final state hadrons 
(cc — ► ng* — > m(qq)). The dominance of annihilation through gluons is most strikingly reflected 
in the narrow widths of charmonium decays into hadronic channels in a mass region where strong 
decays typically have widths of hundreds of MeV. Since the c and the c quarks only annihilate 
if their mutual distance is less than about l/m c (where m c is the c-quark mass) and since the 
average virtuality of the gluons is of the order of 1 — 2 GeV 2 one may indeed expect perturbative 
QCD to be at work although corrections are presumably substantial (m c is not very large). 
In hard exclusive reactions higher Fock-state contributions are usually suppressed by inverse 
powers of the hard scale, Q, appearing in the process {Q = 2m c for exclusive charmonium 
decays), as compared to the valence Fock-state contribution. Hence, higher Fock-state contribu- 
tions are expected to be negligible in most cases. Charmonium decays are particularly interesting 
because the valence Fock-state contributions are often suppressed for one or the other reason. 
In such a case higher Fock-state contributions or other peculiar contributions such as power 
corrections or small components of the hadronic wave functions may become important. Exam- 
ples of charmonium decays with suppressed valence Fock-state contributions are the hadronic 
helicity non-conserving processes like J/ip — > pn or r/ c — ► BB, the radiative J ftp decays into light 
pseudoscalar mesons and , the main topic of this report, the decays of the P-wave charmonium 
states, Xcj- 
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2 Velocity scaling 



Recently, the importance of higher Fock states in understanding the production and the inclusive 
decays of charmonium has been pointed out As has been shown the long-distance matrix 
elements can there be organized into a hierarchy according to their scaling with v, the typically 
velocity of the c quark in the charmonium. One may apply the velocity expansion to exclusive 
charmonium decays as well ||. The Fock-state expansions of the charmonia start as 

\J/il>) = 0(1) |cCi( 3 S x )) + 0(b) |cc 8 ( 3 Pj) g) + 0(v 2 ) |cc 8 ( 3 S!) gg) + 0(v 3 ) 

| xcj) = 0(1) |c Cl ( 3 Pj)> + 0(v) |cc 8 ( 3 S!) g) + 0(v 2 ). (1) 

In |U it is argued that the amplitude for a two-body decay of a charmonium state satisfies a 
factorization formula which separates the scale m c from the lower momentum scales. The decay 
amplitude is then expressed as a convolution of a hard scattering amplitude that involves the 
scale m c , an initial state factor (charmonium wave function) that involves scales of order m c v 
and lower, and a final state factor (light hadron wave function) that involves only the scale Aqcd- 
If m c was asymptotically large, the dominant terms in the factorization formula would involve 
the minimal number of partons in the hard scattering. Terms involving additional partons in 
the initial state are suppressed by powers of v while terms involving additional partons in the 
final state are suppressed by powers of Aqcd/^c- 

The relative strength of various contributions to specific decay processes can easily be estimated. 
Combining the fact that a hard scattering amplitude involving a P-wave cc pair is of order v, 
with the Fock-state expansion ([[]) , one finds for the amplitudes of Xcj decays into, say, a pair of 
pseudoscalar mesons the behaviour 

M( Xc j -» PP) ~ aa 2 s v + ba 2 (v V^F 71 ) + 0(v 2 ), (2) 

where a and b are constants and a s s oft comes from the coupling of the Fock-state gluon to the 
hard process. For the decay reaction J ftp — > BB, on the other hand, one has 

M(J/ij -► BB) ~ aa 3 + ba 3 v(vVa s s oft ) + ca 3 v 2 a s s oft + 0(v 3 ). (3) 

Thus, one sees that in the case of the Xcj the colour-octet contributions are not suppressed by 
powers of either v or l/m c as compared to the contributions from the valence Fock states |J. 
Hence, the colour-octet contributions have to be included for a consistent analysis of P-wave 
charmonium decays. The situation is different for J/i/j decays into baryon-antibaryon pairs: 
Higher Fock state contributions first start at 0(v 2 ). Moreover, there is no obvious enhancement 
of the corresponding hard scattering amplitudes, they appear with at least the same power of 
a s as the valence Fock state contributions. Thus, despite of the fact that m c is not very large 
and v not small (v 2 ~ 0.3), it seems reasonable to expect small higher Fock-state contributions 
to the baryonic decays of the J ftp. 

3 The modified perturbative approach 

Recently, a modified perturbative approach has been proposed Q in which transverse degrees of 
freedom as well as Sudakov suppressions are taken into account in contrast to the standard ap- 
proach of Brodsky and Lepage Q. The modified perturbative approach possesses the advantage 
of strongly suppressed end-point regions. In these regions perturbative QCD cannot be applied. 
Moreover, the running of a s and the evolution of the hadronic wave function can be taken into 
account in this approach properly. 

Within the modified perturbative approach an amplitude for a 2S+1 Lj charmonium decay into 
two light hadrons, h and h! ', is written as a convolution with respect to the usual momentum 



fractions, Xi, x\ and the transverse separations scales, b;, b-, canonically conjugated to the 
transverse momenta, of the light hadrons. Structurally, such an amplitude has the form 



M^( 2S+1 Lj -» hti) = f( c \ 2S+1 Lj) J [dx][dx'] J 





&h denotes the transverse configuration space light-cone wave function of a light hadron. f^ c > 
is the charmonium decay constant. Because the annihilation radius is substantially smaller 
than the charmonium radius this is, to a reasonable approximation, the only information on 

* (c) 

the charmonium wave function required. Tjj is the Fourier transform of the hard scattering 
amplitude to be calculated from a set of Feynman graphs relevant to the considered process, t 

* (c) 

represents a set of renormalization scales at which the a s appearing in T H , are to be evaluated. 
The ti are chosen as the maximum scale of either the longitudinal momentum or the inverse 
transverse separation associated with each of the virtual partons. Finally, exp [— S] represents 
the Sudakov factor which takes into account gluonic corrections not accounted for in the QCD 
evolution of the wave functions as well as a renormalization group transformation from the 
factorization scale \ip to the renormalization scales. The gliding factorization scale to be used 
in the evolution of the wave functions is chosen to be \ip = 1/6 where b = max{bi}. The b space 
form of the Sudakov factor has been calculated in next-to-leading-log approximation by Botts 
and Sterman 

As mentioned before, exclusive charmonium decays have been analysed several times before, 
e.g. [Q]. New refined analyses are however necessary for several reasons: in previous analyses 
the standard hard scattering approach has been used with the running of a s and the evolution 
of the wave function ignored. In the case of the Xcj also the colour-octet contributions have 
been disregarded. Another very important disadvantage of some of the previous analyses is the 
use of light hadron distribution amplitudes (DAs), representing wave functions integrated over 
transverse momenta, that are strongly concentrated in the end-point regions. Despite of their 
frequent use, such DAs were always subject to severe criticism. In the case of the pion, they 
lead to a clear contradiction to the recent CLEO data || on the ttj transition form factor, F^. 
As detailed analyses unveiled, the i ? 7I - 7 data require a DA that is close to the asymptotic form 
oc x(l — x) (?], ||. Therefore, these end-point region concentrated DAs should be discarded for 
the pion and perhaps also for other hadrons like the nucleon f|]. 



According to the statements put forward in Sects. 1 and 2, higher Fock-state contributions are 
neglected in this case. 

From the permutation symmetry between the two u quarks and from the requirement that 
the three quarks have to be coupled in an isospin 1/2 state it follows that there is only one 
independent scalar wave function in the case of the nucleon if the 3-component of the orbital 
angular momentum is assumed to be zero. Since SU(3)f is a good symmetry, only mildly broken 
by quark mass effects, one may also assume that the other octet baryons are described by one 
scalar wave function. It is parameterized as 



4 Results for J/ip (f) -> BB 




(5) 



2 In higher Fock-state contributions additional integration variables may appear. 



channel 


PP 




AA 




AA 


£*£* 




174 


113 


117 


62.5 


65.1 


40.8 


Texp |H] 


188 ± 14 


110 ± 15 


117 ±14 


78 ± 18 


96 ±26 


45 ±6 



Table 1: Results for J/ip — > decay widths (in [eV]) taken from [10|. The three-gluon 
contributions are evaluated with m c = 1.5 GeV, Aqcd = 210 MeV and cib 10 = 0.85 GeV . 
The widths for the decuplet baryon channels are given for one charge state. 



in the transverse momentum space. The set of indices 123 refers to the quark configuration 
u + u_d + ; the wave functions for other quark configurations are to be obtained from (H) by 
permutations of the indices. The transverse momentum dependent part Q is parameterized as a 
Gaussian in fc^/scj (* = 1) 2, 3). The transverse size parameter, ob 8 , appearing in that Gaussian, 
as well as the octet-baryon decay constant, fs 8 , are assumed to have the same value for all octet 
baryons. In j|] these two parameters as well as the nucleon DA have been determined from an 
analysis of the nucleon form factors and valence quark structure functions at large momentum 
transfer (ci£ g = 0.75 GeV -1 ; fs s = 6.64 x 10 -3 GeV 2 at a scale of reference fiQ = lGeV). The 
DA has been found to have the simple form 

^123 A*o) = 60X1X2X3 [1 + 3xi] . (6) 

It behaves similar to the asymptotic form, only the position of the maximum is shifted slightly. 
For the hyperon and decuplet baryon DAs suitable generalizations of the nucleon DA are used. 
The J/i/j decay constant /ju (= /^( 3 Si)) is obtained from the electronic J/ip decay width; it 
amounts to 409 MeV. Except in phase space factors, the baryon masses are ignored and Mj/^ 
is replaced by 2m c for consistency since the binding energy is an 0(v 2 ) effect. 
Calculating the hard scattering amplitude from the Feynman graphs for the elementary process 
cc — > 3g* — * 3(qq) and working out the convolution (Q), one obtains the widths for the J/ip 
decays into BB pairs listed and compared to experimental data in Table |l[ As can be seen from 
that table a rather good agreement with the data is obtained. In addition to the three-gluon 
contribution there is also an isospin-violating electromagnetic one generated by the subprocess 
cc — > 7* — > 3(qq). According to |l0| this contribution seems to be small. 

The extension of the perturbative approach to the baryonic decays of the ip' is now a straight- 
forward matter. One simply has to rescale the J/ip widths by the ratio of the corresponding 
electronic widths 

r( , ^ m = p ^} mB 'y \^f- e+ l\ rm - BB) , (7) 

V ' p p . s .(m B /M J/V) )r(J/^^e+e-) >' KJ 

where p p . s .( z ) = V^- — 4z 2 is the phase space factor. In contrast to previous calculations the 
tp' and the J ftp widths do not scale as (Mj/^/M^,/) 8 since the hard scattering amplitude is 
evaluated with 2m c instead with the bound-state mass. Results for the baryonic decay widths 
of the tp' are presented in Table |2[ Again good agreement between theory and experiment is 
found with the exception of the channel. 



5 Results for Xcj -> PP, J = 0,2 

Let me first discuss the color-singlet contribution to the tttt channel since this is an essentially 
parameter-free contribution. The pion wave function is parameterized analogue to (|) with 
& T = <I>as = 6x(l — x). The pion decay constant, /„•, is 131 MeV and the transverse size 
parameter is given by a n = l/(v / 87r/ 7r which automatically satisfies the ttq — > 77 contraint 



channel 


PP 




AA 




AA 


£*E* 




76.8 


55.0 


54.6 


33.9 


32.1 


24.4 


Texp [ 

[ 


12 
11 


] 
1 


76 ± 14 
53 ± 15 


26 ± 14 


58 ± 12 


23 ±9 


25 ±8 


16 ±8 



Table 2: The three-gluon contributions to the ip' — > decay widths (in [eV]) |l(J. The widths 
for the decuplet baryon channels are given for one charge state. 







PDG [11] 


BES |J 


r[XcO - tt+tt'] [keV] 


45.4 


105 ± 47 


64 ±21 


T[Xc2 -» 7T+7T-] [keV] 


3.64 


3.8 ± 2.0 


3.04 ± 0.73 


r[XcO^^ U ^ U ] [keV] 


23.5 


43 ± 18 




r[Xc2^^ u vr u ] [keV] 


1.93 


2.2 ±0.6 





Table 3: Results for the Xcj decay widths into pions 

(/g) = 1.46 x 10- 3 GeV 2 ) m comparison 

with experimental data. The table is taken from B. 



derived in [13|. This form of the pion wave function is strongly contrained by the ttj transition 
form factor []7j, |8| . The information on the Xcj wave function needed is the singlet decay constant 
which is related to the derivative of the non-relativistic Xcj wave function at the origin 

/ (1) ( 3 Po)=^^=^=^^(0) (8) 
a/3 Vl67rm c 

In -Rp(O) is chosen to be 0.22 GeV 5 / 2 and, as before, m c = 1.5 GeV. The hard scattering 
amplitude is now computed from the Feynman graphs for the elementary process cc — > 2g* — * 
2(qq) and convoluted with the hadronic wave functions according to (||). Evaluating the decays 
width from the singlet amplitudes, one finds, not surprisingly, values which are too small as 
compared to experiment |11]: r(x c o(2) ~~ *■ vr + 7r _ ) = 8.22 (0.41) keV. Of course, the colour-octet 



contribution (see (g)) is lacking and one has to estimate it as well and to add it coherently to 
the colour-singlet amplitude. 

For the calculation of the colour-octet amplitude an octet wave function is required. This 
introduces new degrees of freedom and, for a first estimate, it is advisable to simplify as much as 
possible in order to keep the number of new free parameters as small as possible. It is important 
to realize that in the |ccsg) Fock state not only the cc pair is in a colour-octet state, but also the 
three particles, c, c and g, are in an S'-state. Hence orbital angular momenta are not involved and 
the transverse degrees of freedom can be integrated over. Therefore, one only has to operate 

(S) (8) 

with a DA $V yz\, Z2, 23) and an octet decay constant fj for each J. It seems reasonable 
in the spirit of the non-relativistic expansion to approximate the ccg DA by 6(zs — z) where 
z ~ e/(2m c ) — v 2 /2. e is the binding energy of the charmonium and Z3 refers to the momentum 
fraction the gluon carries. Analogously to the colour-singlet case, the only information on the 

(8) 

colour-octet wave function that then enters the final results is the octet decay constant 
(= /q 8 ^ = . A fit of the unknown decay constant to the experimental data on the decay 
widths provides the value = 1.46 x 10 -3 GeV 2 which appears to be reasonably large |J. 
The resulting decay widths are presented in Table |3[ With regard to the large experimental 
uncertainties fair agreement with experiment is obtained. This analysis can be easily extended 
to other PP channels, see @. 



6 Conclusions 



Exclusive charmonium decays constitute an interesting laboratory for investigating power cor- 
restions and higher Fock-state contributions. In particular one can show that in the decays of 
the Xcj the contributions from the next-higher charmonium Fock state, ccg, are not suppressed 
by powers of m c or v as compared to the cc Fock state and therefore have to be included for 
a consistent analysis of these decays. For J/tp (ip 1 ) decays into BB pairs the situation is dif- 
ferent: Higher Fock-state contributions are suppressed by powers of l/m c and v. Indeed, as 
an explicit analysis reveals, with plausible baryon wave functions a reasonable description of 
the baryonic J/ip (if)') decay widths can be obtained alone from the cc Fock state while for Xcj 
decays into pairs of pseudoscalar mesons agreement with experiment can only be obtained when 
the color-octet contribution is taken into account. 
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